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I. INTRODUCTION 

Charge transport through topologically nontrivial ma- 
terials is of great current interest, offering novel fun- 
damental insights as well as potential applications in 
topological quantum computing.— The recently discov- 
ered topological insulators as well as topological super- 
conductors (TSs^PEl are predicted to exhibit spectacu- 
lar nonlocal transport phenomena and have consequently 
attracted a lot of attention. In particular, for one- 
dimensional TS quantum wires, the crucial role of Ma- 
jorana bound states (MBSs) located near the interface 
to topologically trivial regions has been emphasized 
Majorana fermions are special in that they are their own 
antiparticles, i.e., the fermion creation operator is equal 
to the annihilation operator, and the condensed-matter 
setup discussed below could offer experimental signatures 
for these elusive and hitherto unobserved particles. 

We here study nonlinear electronic transport through 
a mesoscopic TS wire (also referred to as "dot") contain- 
ing a pair of MBSs. The wire is assumed to be con- 
tacted by conventional normal-conducting leads. The 
contacts correspond to thin insulating layers or barriers 
and will be modeled by a tunnel Hamiltonian. We include 
Coulomb interactions on the dot via the charging energy 
E c = e 2 /(2C). The schematic setup studied in this work 
is shown in Fig. [T] and we briefly outline possible exper- 
imental realizations below. The noninteracting (E c = 0) 
version of this setup was studied in detail before EEHIU 
and one finds resonant Andreev reflection with the max- 
imal (unitary) value Gl/h = 2e 2 /h of the linear conduc- 
tance through the left/right contact. For the interacting 
case {E c > 0), an important work by FrP^ has provided 
a general framework to address transport in this setup. 
The Hilbert space constraint introduced in Ref. 13 is, 
however, difficult to handle in actual calculations. In fact, 
the concrete results presented in Ref. |T3] (see also below) 
were mostly restricted to the strong Coulomb blockade 
(CB) regime realized for large E c . We formulate an alter- 
native Keldysh functional integral approach below, where 
the constraint is automatically satisfied. We give detailed 
results for the weak CB (small E c ) regime, where we find 
a pronounced interaction (E c ^ 0) induced suppression 




Figure 1: (Color online) Setup studied in this work: a topo- 
logical superconductor (TS) dot containing Majorana bound 
states at the ends (indicated by black dots) is coupled via tun- 
nel contacts (hybridizations Vl,r) to normal-conducting (N) 
metallic electrodes. The charging energy E c = e 2 /(2C) causes 
phase fluctuations on the TS dot and influences the current- 
voltage relation. We assume that the superconducting gap A 
of the TS is large compared to the other energy scales, i.e., 
quasiparticle excitations on the dot can be neglected. 



of the current. 

Experimental realizations for the setup shown in Fig.[T] 
arise under several distinct physical scenarios. For in- 
stance, a nanowire made of a 3D topological insulator 
(e.g., Bi 2 Se 3 ) with proximity-induced superconducting 
correlations is predicted to have MBSs at both endsfHl 
Similar setups can be realized in 2D topological insu- 
lators with "quantum spin Hall" edge states! 13 * 15 ^ An- 
other example is based on the vortex core states in 
p-wave superconductors^ or in 2D noncentrosymmetric 
superconductors.^ Finally, the MBS should also be re- 
alizable in semiconductor quantum wires with a Zee- 
man field and proximity-induced superconductivity.^ 
The stability of the MBS when including interactions in 
an isolated wire has also been addressed in several recent 
works. — However, so far no experimental signature for 
the MBS has been reported in any system. 

The structure of this paper is as follows. In Sec. [n] 
we describe the model employed in our study and briefly 
show that our formulation reproduces known results^ 
in the strong CB limit. We then construct a general 
Keldysh functional integral representation for the inter- 
acting problem in Sec. |III| This theory is employed in 
Sec. |IV| to evaluate the current- voltage characteristics un- 
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der weak CB conditions. We briefly summarize our main 
findings and provide an outlook in Sec. ^ Technical de- 
tails and derivations can be found in three appendices. 
Note that we often use units where h = kg = 1. 



II. MODEL 



fermion parity of the dot is determined by rid only and 
not affected by changes of N. By construction, both N 
and x commute with jj, and therefore no constraint on 
the Hilbert space arises. 

The full Hamiltonian describing transport through the 
dot, H = H c + H t + Hi, contains the Coulomb charging 
term 



We consider transport through a mesoscopic TS quan- 
tum wire contacted by normal-conducting leads, see 
Fig. [T] We assume throughout this work that the su- 
perconducting gap A is the largest relevant energy scale, 
and thus no quasiparticles above the gap are involved. In 
particular, max(eV, E c ) < A. We study a "floating" (not 
grounded) TS dot, where current must be conserved un- 
der steady-state conditions. In some previous works, see, 
e.g., Ref. |9j a grounded superconductor was considered, 
where current need not be conserved. Note also that 
the situation is different for conventional (topologically 
trivial) superconducting dots, where transport is typi- 
cally mediated by cotunneling processes and the resulting 
conductance is much smaller.^ In the present case, two 
MBSs may host a fermion at no energy cost, and thus no 
even-odd parity effeclP^ arises. 



A. Model 

For realistical TS wire length, the two MBSs located 
near the left/right (j = L/R) contact correspond to de- 
coupled Majorana fermion operators, jj — 7J, with an- 
ticommutator relation {7i,7j} = <%• In the low-energy 
sector of interest here, only the MBSs provide available 
fermion states to tunnel through the superconductor. It 
is convenient to introduce the complex auxiliary fermion 
d = {il + ryflVv^, i-e., 

lL = (d + S)/V2, lR = -t(d-d^)/V2. (1) 

For an isolated wire, the fermion parity (— ) n of the two 
degenerate ground states with even (n = 2N) or odd 
(n = 2N + 1) number of electrons fixes the occupation of 
the d fermion level, 



H r 



E C (2N + n d - n ) 2 



(3) 



2J7L7-R = 2n d -l= (-) 



n+l 



rid 



dU. 



(2) 



This imposes a constraint on the Hilbert space, 13 and the 
Majorana fermion dynamics cannot be taken as indepen- 
dent of the charge (or the dual phase) dynamics. Below, 
we instead formulate an alternative but essentially equiv- 
alent approach free of any Hilbert space constraint, see 
also Ref. which is technically easier to handle. To 
that end, consider the number operator N for Cooper 
pairs conjugate to the condensate phase x on the dot, 
[x, N] — i- The dot's instantaneous charge state is then 
fully determined by specifying the configuration (N,n d ), 
where N g No and rid = 0, 1 are eigenvalues of N and 
h d , respectively. Importantly, Eq. ^ tells us that the 



where no G K is tunable by the backgate voltage V g in- 
dicated in Fig. [l] Electrons in the leads correspond to 
effectively spinless fermion operators Cj = L/R^ with mo- 
mentum k, where for given j, Cjk and Cj k are tunnel- 
coupled to the respective Majorana fermion 7, only. For 
simplicity, we assume a /c-independent tunnel matrix el- 
ement Xj encoding the overlap of the lead state with the 
respective MBS wavefunction. As detailed in App. [A] 
the tunnel Hamiltonian then takes the form 



- iX 



fl4 fc (d-e-*dt) 



~ lx d)) 
+ h.c, 



(4) 



which is charge conserving. In particular, the differ- 
ent terms in H t describe electron tunneling from the 
dot to lead j, either by destroying the d state without 
changing N, i.e., (N, 1) — > (N,0), or by occupying the 
d state together with a splitting of a Cooper pair, i.e., 
(N, 0) — > (N — 1,1), plus the conjugate processes. As 
elaborated in Appendix[B] we treat the lead Hamilt onian 
Hi within the standard wide-band approximationJ^ilsU 
The applied bias voltage V corresponds to the chemical 
potential difference in both leads, eV = (jll — n R- For 
further reference, we also define the hybridization energy 
scales 



r, = 2™, I a, | 2 , 



(5) 



where Vj is the density of states in lead j = L/R. 

It is now straightforward to rederive the well-known 
results®021 f or t ne noninteracting (E c = 0) limit from the 
above Hamiltonian, see also Sec. |IV| Before turning to 
the construction of the Keldysh functional integral for 
the interacting problem in Sec. |III| let us briefly show 
that the above formulation easily reproduces previous 
results 13 reporting electron "teleportation" under strong 
CB conditions, i.e., for max(eV, Tl,r) <C E c . 



B. Strong Coulomb Blockade 

Under strong CB conditions, the charging energy ^ 
is dominant and needs to be examined first. We focus on 
the resonant situation realized for half-integer values for 
no, where two different cases arise: (i) For no = 2£+ 1/2 
with integer I, degeneracy is achieved for fixed Cooper 
pair number N = £, where the d state is either occupied 
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(rid = 1) or empty (rid = 0). (ii) When n = 21 - 1/2, 
the degeneracy point is achieved for different Cooper pair 
numbers, namely for (N, n^) = (£ — 1, 1) and (£,Q). We 
conclude that the low-energy description needs to keep 
only two states, either of type (i) or (ii) depending on 
no. For case (i), N remains constant and all terms cx e ±tx 
(where N is raised or lowered by one unit) can be omitted 
in Eq. Q. The tunnel Hamiltonian then takes the form 

H t ] = 4 E (> C L iAflC^) d + h.c. 

h 

Under scenario (ii) only transitions (I — 1,1) <H- (£, 0) 
are possible. Using / = e~ %x S as effective single-charge 
fermion, we arrive at 

H < U) = 71 E ( A * C L + iA «4, fc ) f + t-c 

In both cases, we recover the resonant tunneling model 
describing teleportationP^ This leads to the unitary con- 
ductance value Gl = Gr = e 2 /h instead of the resonant 
Andreev reflection value 2e 2 /h found when E c = 0. 

III. KELDYSH FUNCTIONAL INTEGRAL 

We now turn to the general Keldysh functional 
integrapil22l f or the interacting problem. This corre- 
sponds to the Majorana generalization of the (real-time 
version of the) Ambegaokar-Eckern-Schon (AES) action 
for a conventional metallic grain pHEU 

A. AES action for Majorana induced transport 

The relevant low-energy degree of freedom character- 
izing the interacting problem is the single-electron (half 
Cooper pair) counting phase field cj>(t) = %/2. For the 



nonequilibrium problem at hand, we employ the text- 
book Keldysh functional integral formulation^ 21 * 22 ! and 
double the phase field <f)(t) — > 4>±(t) according to the 
forward and backward branch of the Keldysh time con- 
tour. Technically, the phase field is introduced through 
a Hubbard-Stratonovich transformation of the charging 
energy in Eq. ([3]), see App.[B]for details. Physically, <fi is 
the dual variable to the Cooper pair charge N. We then 
introduce classical/quantum phase fields 1241251 

= ^ (>+ + </>-), 0,(t) = 0+-0- (6) 
The charging energy ([3| yields the action contribution 

S c = Jdt^ q (^-+n \ (7) 

encoding Coulomb blockade effects. By virtue of the 
Hubbard-Stratonovich transformation, the fields corre- 
sponding to the lead fermions and to the Majorana 
fermions effectively become noninteracting. In the func- 
tional integral, those fields can therefore be integrated 
out analytically; for technical details, see Appendix [5] 
The final (phase representation of the) functional inte- 
gral for the Keldysh partition function reads 

Z = J V(<j> c ,<j> q ) e l ^ +s '\ (8) 

The fermion field integrations result in the action piece 

S f = -i lnPf [rxiW-O+^-feO] ' ( 9 ) 

j = L/R 

where "Pf" denotes the Pfaffian and f x ,y, z are Pauli ma- 
trices in (rotated^!) Keldysh space. Moreover, we define 
the antisymmetric (time and Keldysh) matrix 



Aj(t,t') = 2iY j F(t-t') 



c(t)c(f) cos $j (t, t') c(t)s(t') sin $j (t, t') 
-s(t)c(U) sin t') s(t)s(t') cos $,•(*, f) 

c 2 (t)5-{t-t')-s 2 (t)5+{t-t') 



(10) 



-c 2 (t)5 + (t-t') + s 2 (t)5_(t-t') 







with c(t) = cos[<p q (t)/2],s(t) = sm[(/) q (t)/2], and 8±(t) = 
5(t±0 + ), where the infinitesimal shifts reflect the proper 
causality features. At temperature T, the distribution 
function F(t — t') has the Fourier transform 



F(e)=tanh[e/(2T)], 



(11) 



and the phase function appearing in Eq. ( 10 1 is 

^(MO = - + 4>c(t) - 4>o(t'). (12) 

Finally, the hybridizations Tj were defined in Eq. (J5|. 

While Eq. with the action in Eqs. Q and (JsJ , pro- 
vides a general representation of the interacting nonequi- 
librium Majorana transport problem, approximations are 
necessary in order to obtain concrete analytical results. 
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B. Semiclassical expansion 

We here are mostly interested in the weak CB regime, 
where standard arguments^ imply that a semiclassical 
approximation in the phase representation is appropriate. 
In particular, fluctuations of <p„ around zero are small, 



i.e., we can expand the action in Eq. ([. 
order in <p q . Keeping the full nonlinear 



to quadratic 
, dependence, 

the matrix ( 10 1 is thus approximated by A J=i / R ~ A^ -* + 

' (t, t') of order 



, Aj 2 " 1 . Specifically, the matrices A^ r ' 



A (o) 

A 
A 



j 

(i) 

j 

(2) 



( 2F(t-i?)ca8$ j (t,tf) 6-(t-tf) 



-5+(t-t>) 



Y f ^2 

-i-^F(t - t') cos $j(i,t') 



-2^(t)0,(f) 



±2 



K{t)5{t-t') 



At this stage, it is convenient to introduce the interacting 
Keldysh-Majorana Green's function (GF) 



GJt,t') 



f x id t + A 



(o- 



(13) 



This has a triangular representation in (rotated) Keldysh 
/ o G A \ 

space, Gj — ( nR „^ I , with retarded/advanced GFs 

G^ A and Keldysh GF G^ . Note that in this representa- 
tion, the Majorana GF has a bosonic Keldysh structure, 
see also Ref. IT21 According to Eq. (13 1, G^ A obeys the 



noninteracting (ip c = 0) equation of motion and thus has 

the Fourier representation Gj (e) = l/(e± iTj). This 
implies the spectral function 



A i (e) = -Im[Gf(e)] = 



e 2 + rj 



(14) 



Only G^ is affected by interactions encoded in cj) c . 
The action ^ now takes the form 

S f = -iJ2 lnPf {Gj 1 ) + Sf ] + Sf + 0(<f> 3 q ), (15) 



where the first (<f> q — 0) term vanishes as a result of the 
unitary time evolution along the closed Keldysh contour. 
The first-order (in <f> q ) contribution is 

S ? = - 2 E Tr ( 6 i A f = / dt Wo ( 16 ) 



where the trace "Tr" extends both over time and Keldysh 
space, and 1 denotes the total current flowing into the 
dot. Some algebra yields 



X(t) = t j J ^ Gf{t - t') {F jt .(lf - t) cos [<t> c (t') ~ M*)} + p jAt' ~ *) sin [<£ c (t') - 4> c {t)]} (17) 



with the time-symmetric (-antisymmetric) functions 



IV. WEAK COULOMB BLOCKADE LIMIT 



Fj >s {t) = F(t) wn(njt), F Ao (t) = F(t) cosQijt). (18) 



A. Langevin approach 



Finally, the second-order term in Eq. ( 15 ) reads 



(19) 



When the charging energy E c is sufficiently small, fluc- 
tuations of the phase fields <j> c ^ q are small and allow for 
a quadratic expansion of the action Sf in both fields.^ 
This expansion is detailed in Appendix [C] and allows to 
reformulate the functional integral (|8| in terms of the 



5 



The damping kernel rj(t — t') has the Fourier representa- 
tion 

(20) 



f,( W )= Yl 77 / ^(Ajie + Ljm+Ajie-LjMiF^e + uM-Fj^e-wp)]. (21) 

The Gaussian noise field has zero mean value and the correlation function (£(i)£(t')) = — t'), where the 
fluctuation kernel has the Fourier representation 

K(oj) = Ey/ ^(^(e + ^/2) + ^.(e-o;/2))[l-^ )0 (e + W /2)F j)0 (e-a;/2^ (22) 

- E r ' / ^ Re + «V 2 ) G ?( e -w/2)]F ijS ( e + «/2)^,.(c -«/2). 

I 

In equilibrium {hl,r — 0), Fj s = and Fj a — F, and 
hence the fluctuation-dissipation relation 

K cq (uj) = | coth ?7 cq (a;) (23) 

is satisfied, providing an important consistency check. 

In the zero-temperature limit, both rj(u>) and K ((d) can 
be evaluated in closed form. Specifically, we find for the 



equivalent semiclassical Langevin equation, 

^Mt)+ f dt' »7(t-O0c(O = C(*)- 



damping kernel 




while the fluctuation kernel is 



K T=0 (u) 



\u\t)(uj) 



E 



^e(2| Mj | 



M) 



tan 



r,- 



tan 



r i]n (N-M) 2 + r 



(25) 



where O is the Heaviside function. 



The Langevin equation ( 20 ) allows for an intuitive in- 



terpretation of the charge dynamics in this interacting 
Majorana system in terms of an effective RC circuit. Its 
solution is used below in order to compute the current- 
voltage relation. 



B. Current-voltage relation 

Next we outline the calculation of the dc current. 



I = II = —Ir. By using Eq. (17) and the noise av- 
erage (■■•)? — / (' ' ' )e~^ K ? j we find the (indeed 
t-independent) result 



dt 1 Gf(t - t') F hS {t' - i)(cos[0 c (i') - c (i)]) 4 



+ Fj >a (t' - i)(sh# c (i') - Mt)])i 



Fj,a/s is given in Eq. (18 1, and 



4>c(t)= j'dT D R {t-T)i(T) 



(26) 



solves Eq. (20 1 for a given noise trajectory, where D R (t — 
t) is the retarded solution of Eq. (20 1 with — > S(t— r). 
Performing the Gaussian noise average, we obtain 



dt' Gf(t - t')F jtS (t' - t)e- J{t ~ l '\ (27) 



where the noise correlations are encoded in J(t) = J(—t). 
similar to the standard U P(E)" theory of dynamical 

CB m 

J(t-t')= l -([Ut)-Ut')] 2 ) i - (28) 

Note that in order to obtain J(t—t') one first has to solve 
the Langevin equation ( 20j) for a given noise trajectory, 
and then average over all noise realizations. 
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Since J(t — t') > 0, interactions always decrease the 
current (27 1 with respect to the noninteracting solution. 



lf\ The latter follows by setting J -> in Eq. ((27 1 



In fa ct, after Fourier transformation, we get the known 
resultjUSm] 



/. (0) = r, 



de 

2n 



(29) 



Since 770 < 4/7T, see Eq. ([Ml, and using E c = e 2 /(2C), 
the effective resistance R in Eq. ( 34 I always fulfills R > 
h/(4e 2 ). Note that this corresponds to a parallel resis- 
tance as seen by the dot. 

However, at the same time it is essential to keep the 
full frequency-dependence in the fluctuation kernel K(w) 
when evaluating J(t) in Eq. ( 28 1 . Using Eq. ( 33 1 , some 
algebra yields 



with the Lorentzian spectral function in Eq. ( 14 1 . In lin- 
ear response (/x, — > 0) and taking the zero-temperature 
limit, we recover resonant Andreev reflection with quan- 
tized conductance Gj — edlj/d/ij — 2e 2 /h. 

Equation (27) for the current- voltage relation in the 



weak Coulomb blockade regime is one of the main results 
of this work. We use this expression later on to quan- 
tify the effect of Coulomb interactions (E c =^ 0) on the 
transport properties in this Majorana system. 



Retardation effects 



The Langevin equation ( 20 ) in general contains mem 



ory effects because of the frequency-dependence of the 
damping kernel. These are most pronounced for T = 0, 
where r)(uj) takes the form in Eq. (24). For u — > 0, this 
gives the finite value 



J(i) 



nVo Jo 



w 2 (l + w 2 /fi 2 )' 



(35) 



where the current I = II follows from Eq. (27 1 



D. Zero-bias anomaly 

We now present results for the IV relation in the 
zero-temperature limit, where MBS effects are most pro- 
nounced. We adopt the Markovian approximation de- 
scribed in Sec. |IV Cj and study a symmetric system, where 
= Tr = r/2 and [i^ — — /Ur — eV/2. As relevant 
transport quantity, we define the dimcnsionlcss nonlin- 
ear conductance 



g{V) 



i(v) 

e 2 V/h' 



(36) 



1 



m = 



TT ^ 1 

j 



(30) 



t'l 



while for — > 00, we have r)(uj) ~ r/|w| with r = Tl 

Tr. On low frequency scales, |w| < min {^\J 1 

the kernel rj(ui) can be approximated by r]{ui) ~ 770- This 
approximation corresponds to the absence of retardation 
in the damping kernel rj(t). This approximation works 
even better for finite T, where 



Vo 



j 



de 



1 



r. 



T J 2^ cosh 2 (^) (e-Mi) 



(31) 



yields a damping suppression, 770 ~ T/(2T), in the high-T 
limit. 

We therefore employ the simplification rj(uj) = t]q 
from now on, which is appropriate to describe low-energy 
transport. This implies from Eq. (20) that we now have 
a Markovian Langevin equation, 



4> c {t) + Slj> c (t) = 2E c £(t), 



and D R (t) in Eq. p6l) simplifies to 



D R {t) = 1 (1 

m 



-nt\ 



8(t). 



The inverse U RC time" in Eq. (32 1 is given by 
Q = {RC)- 1 = r l0 E c . 



(32) 



(33) 



(34) 



For E c = 0, the exact solution [Eq. (29)] yields 



g(0) {V) = J_ tan -l'' r 

ev 



(37) 



Note that the unitary limit (resonant Andreev reflection) 
corresponds to g = 1, which is reached for V — > in 
Eq. §7}. 

We now write J(t) = J + J 1; where 



1 

27T7/0 



duj 1 — cos(cjf) 

V 1 + cu 2 /n 2 



(38) 



comes from the first term (oc rj(u>) = ?7o) in Eq. (25 1. Ji is 
a pure nonequilibrium term: the relevant kernel contribu- 
tion scales as ~ V 3 . Let us then first consider the linear 
response regime, where J\ can be discarded. To logarith- 
mic accuracy, Jo(t) ~ (27r?7o) _1 rn(f2t) for fit > 1, where 
770 = 4/7r within linear response, see Eq. (30 1. Equation 



( 27 1 then yields the dimensionless linear conductance in 



analytical form, 



g(V -»• 0) = coiEjT)- 1 / 8 , c » 0.96. 



(39) 



This result hokLpS for V < E c . As a function of E c /T, 
Eq. ( 39 ) reveals a power-law suppression of the linear 
conductance with the universal exponent 1/8, reminis- 
cent of a zero-bias anomaly. For arbitrary E c /T, the cur- 
rent given by Eq. ( 27 1 can be evaluated numerically. The 
resulting linear conductance is shown in Fig. [2] and the 
analytical result ( 39 1 is quite accurate in describing the 
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Figure 2: (Color online) Dimensionless linear conductance, 
9 — g{V ~ ► 0)j vs E c /T in the weak CB regime. The black 
dashed curve has been obtained from numerical analysis of 
Eq. (p7|, with J(t) in Eq. |35| and K(u) in E q. |25| . The 



red solid curve gives the analytical result in Eq. ( 39 1 valid for 
E c > r. 



E c /T=0.157 

eV=1.57 

E°=0 




0.4 



ev/r 

Figure 3: (Color online) Dimensionless nonlinear conduc- 
tance g(V) vs eV/T for several values of the charging energy. 
Results for E c > were obtained by numerical integration, 
the result for E c — is exact. 



interaction (E c ^ 0) induced suppression of the current 
for E c > T. 



Next we turn to the nonlinear conductance g(V) in 
(36 1. Numerical results for several values of E c /T 



lown in Fig. [3] where the exact result [Eq. ( 3T l] for 
is also depicted. We find a clear suppression of 



Eq 
are s 
E c = 

g(V) when the applied bias voltage increases. This sup- 
pression becomes stronger with increasing E c . In fact, we 
observe that the interaction induced suppression of the 
current is more pronounced for small voltage. We there- 
fore interpret our results as a zero-bias anomaly caused 
by the (weak) Coulomb blockade of Majorana fermion 
induced transport. 



V. CONCLUDING REMARKS 

In this work, we have presented a general theory of 
transport through a mesoscopic superconductor contain- 
ing a pair of MBSs. We have studied the regime of low- 
energy transport, where metallic leads are attached and 
all relevant energy scales are small compared to the su- 
perconducting gap. The case of weak Coulomb blockade, 
where the charging energy E c does not significantly ex- 
ceed the typical hybridization scale T, has been investi- 
gated in some detail, and we found clear signatures for 
an interaction (E c ^ 0) induced suppression of both the 
linear and the nonlinear conductance. 

We have introduced a formulation to analyze interac- 
tion effects in electronic transport through MBSs which 
does not rely on the Hilbert space constraints of previ- 
ous approaches and thus is easier to handle in general. 
We are presently extending this formulation to study 
several other situations of current interest. In particu- 
lar, relevant questions include the full counting statis- 
tics and the nonequilibrium dephasing in multi-terminal 
transport geometries containing MBSs, and the Joseph- 
son effect with superconducting leads. Furthermore, the 
strong Coulomb blockade regime also allows for the de- 
tailed study of the full Keldysh functional integral (|8| in 
terms of a rate equation approach. 

To conclude, we are confident that the characteristic 
features reported here, i.e., the suppression of the lin- 
ear conductance with increasing E c /T and of the nonlin- 
ear conductance with V/T in the limit of weak Coulomb 
blockade, can be observed experimentally once MBSs 
have been detected. The setup studied here may al- 
low to find signatures of the elusive Majorana fermions 
in transport properties, including phenomena such as 
resonant Andreev reflection, teleportation, and the uni- 
versal power-law suppression of the linear conductance, 
g ~ (£ , c /r) -1 / 8 , for intermediate values of the ratio 
E c /T. 
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Appendix A: Derivation of the tunnel Hamiltonian 

In order to keep the paper self-contained, we here pro- 
vide a derivation of the tunnel Hamiltonian Q. Follow- 
ing FlensbergJS] the field operator ^a-(x) for electrons 
with spin projection a ^ in a TS wire of length L 
(with < x < L) defines the Majorana fermion opera- 
tors 



lj=L/R 



E 
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with the MBS wavefunction fj a (x). Projecting the full 
Nambu spinor onto the MBS subspace, one obtains^ 

3 

The tunnel Hamiltonian now reads 




where Ct fe(j creates an electron in lead j = L/R with 
momentum k and spin projection <r. Here tj(x) is as- 
sumed fc-independent for simplicity, which defines the 
tunnel matrix elements Vj a = J dx tj(x)f* (y (x). Note 
that one can always form suitable linear combinations of 
Cjk,-f and Cjk,i to form spinless lead fermions Cjk coupled 
to the MBSs,' 

The other (orthogonal) linear combination then decou- 
ples from the problem. Finally, inserting the auxiliary 
fermion representation (|T|) and taking into account charge 
conservation, we arrive at the tunnel Hamiltonian quoted 
in Eq. @. 

Appendix B: On the Keldysh functional integral 

Here we provide a detailed derivation of the AES action 
in Sec. |III A| We start by constructing the Lagrangian 
L c for the isolated dot in terms of <f>(t) and the Grass- 
mann variable d(t), plus the corresponding "velocities" 
4> = \dH c /dN and d. Noting that for Grassmann vari- 
ables, id is canonically dual to d, we have 

ft ■ - 

L c = — + n 4> + d(id t - 4>)d- 

Adding the tunnel contribution and performing the gauge 
transformation d — > e~ ll ^d, we obtain 

<t> 2 ; iv^ ■ 

Lc + L t = +n a (f>+ -^Ijlj 

c j 

E( A w^s+h.c). 



Here we switched back from the auxiliary d fermion to the 
Majorana field jj. Note that up to a full time derivative, 
idd —> S,-7j7j- The Grassmann fields {ipjk,^Pjk) 

correspond to the lead fermion operators (c . fc , ct fc ). Using 
the Keldysh formulation, we double all fields according 
to the forward and backward branch of the Keldysh time 
contour, i.e., jj(t) — > ("fj,+ (t),jj,-(t)) T and so on. It is 
also convenient to gauge out the chemical potentials Lij 
in the leads, ipjk,±{t) —> e'^'^t^ (t). We now use the 
Keldysh matrix notation (p(t) — 4> c + f z (/) q /2, see Eq. pj. 
The complete Keldysh action, S = SL + S c + S t + Si, 
contains S c in Eq. ^ and the pieces 

S ~< = lYlj dt H f *% 

S * = _ £ A */ ^^e- i ^ t +^7 J -+h.c, 

jk 

Si = ^2 J ' dt ipjkf z {idt- €jk)i>jk, 

jk 

where tjk refers to the dispersion relation in lead j. 

The current Ij flowing from lead j into the dot fol- 
lows from the Heisenberg equation of motion. The cur- 
rents obey the relation (ijr + Iii)(t) = SSt /5<f> q (t). Cur- 
rent conservation implies (Ii) = —{Ir), which fixes the 
mean value (ft) to the chemical potential ii s of the su- 
perconducting dot. The latter has to be determined self- 
consistently from current conservation. We then redefine 
(j) c — > Li s t + 4> c {t), i.e., <fic(t) now refers to fluctuations 
around the mean-field value (/x s t) of the classical phase 
variable. The ensuing changes in S c and S t can be ab- 
sorbed in a renormalization of no and 11 j. In particular, 
denoting their bare values by fio and flj, respectively, we 
find no = Hq + fi s /(2E c ) and Lij = p,j + fx s . 

The next step is to integrate ou t the Gras smann fields 
(ipjk, 4>jk)- This is a standard step 1 21 ' 22 ' 24 ' 2 ^! and leads to 
the effective Majorana action S c g replacing Si + St + 5 7 . 
Before turning to the result, we perform the usual ro- 
tation in Keldysh space in order to have triangular 
Green's function (GF) representations, using the unitary 

matrix-^!/ = ( J ^ \ . The rotated Majorana fields 

are 7,- = Ljj. Evaluating the resulting momentum inte- 
grals for the leads in wide-band approximation, we find 
the effective Majorana action 
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-<5 + (t-t') 



The hybridizations are defined in Eq. (|5| with lead 
density of states Va = J2k^( e jk), F(t) has the Fourier 



Eq. (12) 



transform (111, and the phase function $ , is specified in 



Now we are ready to also integrate out the Majorana 
fields. This finally yields Eq. ^ where 



where r]{t — t') is determined by Eq. (21 1. Next we note 
that S^p , which is already of order <fi q , can be evaluated 

by replacing the interacting Keldysh GF with its 
noninteracting (d> c = 0) version. After some algebra, we 
find from Eq. <[l9l 



is the antisymmetric part of Qj. 



Appendix C: Derivation of Langevin equation 

Here we provide some details concerning the deriva- 
tion of the Langevin equation p0| in Sec. IV A We first 

analyze the action piece Sf — j dt T<j) q . As mentioned 
in App. [B] the chemical potential /i s of the dot has to 
be chosen self-consistently in order to fulfill current con- 
servation. By construction, we must therefore have the 
average value 1 = when 4> c — 0. Using Eq. (17 1, ex- 
panding to first order in <j) c . and performing a partial 
integration, we find 



dt <f> q (t) / dt' r)(t-t')<p c (t'), 



Sf = l - / dtdt' <f> q (t)K(t - t')<t> q {t') 



with K(t - t') determined by Eq. (|22j) . The 
Fourier-transformed kernel K(ui) directly describes finite- 
frequency current noise correlations for E c — 0, a ques- 
tion that has been studied for to — > in Ref. [TU Finally, 
we complete the derivation by Hubbard-Stratonovich 
transformation of the <f> 2 q term in the generating func- 
tional [Eq. M], 



e iS f 



I 



PC e 



-i/dtdt' £(i)K _1 (t-t')S(t')+iJ<i* 4><,t 



Functional integration over </>„, which appears only lin 



early in the action, yields Eq. (20 1, and the Gaus 



sian noise field Hi) indeed has the correlation function 
K(t-t') in Eq. p2l). 



1 C. Nayak, S.H. Simon, A. Stern, M. Freedman, and S. Das 
Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

2 M.Z. Hasan and C.L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

3 X.L. Qi and S.C. Zhang, Rev. Mod. Phys. (in press); 
larXiv: 1008.20261 

4 A.Yu. Kitaev, Phys. Usp. 44, 131 (2001). 

5 L. Fu and C.L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

6 Y. Tanaka, T. Yokoyama, and N. Nagaosa, Phys. Rev. 
Lett. 103, 107002 (2009). 

7 K.T. Law, PA. Lee, and T.K. Ng, Phys. Rev. Lett. 103, 
237001 (2009). 

8 J. Linder, Y. Tanaka, T. Yokoyama, A. Sudt>0, and N. 
Nagaosa, Phys. Rev. Lett. 104, 067001 (2010). 

9 C.J. Bolech and E. Dernier, Phys. Rev. Lett. 98, 237002 
(2007). 

10 J. Nilsson, AR. Akhmerov, and C.W.J. Beenakker, Phys. 
Rev. Lett. 101, 120403 (2008). 

11 K. Flensberg, Phys. Rev. B 82, 180516(R) (2010). 

12 A. Golub and B. Horovitz, Phys. Rev. B 83, 153415 (2011). 



13 L. Fu, Phys. Rev. Lett. 104, 056402 (2010). 

14 M. Cook and M. Franz, |arXiv:1105.1787| 

15 L. Fu and C.L. Kane, Phys. Rev. B 79, 161408(R) (2009). 

16 M. Sato and S. Fujimoto, Phys. Rev. B 79, 094504 (2009). 

17 R.M. Lutchyn, J.D. Sau, and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010); Y. Oreg, G. Refael, and F.v. 
Oppen, ibid. 105, 177002 (2010). 

18 S. Gangadharaiah, B. Braunecker, P. Simon, and D. Loss, 
Phys. Rev. Lett. 107, 036801 (2011); E.M. Stoudenmire, 
J. Alicea, OA. Starykh, and M.P.A. Fisher, Phys. Rev. 
B 84, 014503 (2011); E. Sela, A. Altland, and A. Rosch, 
Phys. Rev. B 84, 085114 (2011). 

19 D.V. Averin and Yu.V. Nazarov, Phys. Rev. Lett. 69, 1993 
(1992). 

20 B. van Heck, F. Hassler, AR. Akhmerov, and C.W.J. 
Beenakker, |arXiv:1108.T095l 

21 Yu.V. Nazarov and Ya.M. Blanter, Quantum transport 
(Cambridge University Press, Cambridge, 2009). 

22 A. Altland and B.D. Simons, Condensed matter field the- 
ory, 2nd edition (Cambridge University Press, Cambridge, 



10 



2010). 

23 G. Schon and A.D. Zaikin, Phys. Rep. 198, 237 (1990). 

24 Yu.V. Nazarov, Phys. Rev. Lett. 82, 1245 (1999); D.S. 
Golubev and A.D. Zaikin, Phys. Rev. Lett. 86, 4887 
(2001); Phys. Rev. B 69, 075318 (2004); P.W. Brouwer, 
A. Lamacraft, and K. Flensberg, Phys. Rev. B 72, 075316 
(2005). 

25 A. Altland and R. Egger, Phys. Rev. Lett. 102, 026805 



(2009). 

H. Grabert and M.H. Devoret (eds.), Single Charge Tun- 
neling, NATO ASI Series B, Physics v.294 (Plenum Press, 
New York, 1992). 

The weak CB regime is eventually left for P <C E c , where 
Eq. (391 ceases to be valid. 



